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Abstract: Image registration has great significance in medicine, 
with a lot of techniques anticipated in it. This research work 
implies an approach for medical image registration that 
registers images of the mono modalities for CT or MRI images 
using Modified Adaptive Polar Transform (MAPT). The 
performance of the Adaptive Polar Transform (APT) with the 
proposed technique is examined. The results prove that MAPT 
performs better than APT technique. The proposed scheme not 
only reduces the source of errors and also reduces the elapsed 
time for registration of brain images. An analysis is presented 
for the medical image processing on mutual- information-based 
registration. 

Index Terms: Adaptive Polar Transform; Gabor wavelet 
transform; Modified Adaptive Polar Transform; Modified 
Gabor wavelet transform; Registration; Brain Images. 

I.INTRODUCTION 

Image registration is the method of superimposing 
images of the similar view captured at various times, from 
different perspectives, and/or by different sensors. The 
registration geometrically aligns two images (the model and 
sensed images). These approaches are classified based on 
their nature (area based and feature-based) and have four 
fundamental steps: (i) feature detection, (ii) feature matching, 
(iii) mapping function design, and (iv) image transformation 
and sampling. It is essential in all image analysis tasks in 
which the ultimate information is realized from the grouping 
of a range of data sources like in image fusion, change 
detection, and multichannel image restoration. Normally, 
required in medicine (combining CT and / or MRI data to 
obtain additional complete information about the patient, 
monitoring tumour growth, treatment verification, evaluation 
of the patient's data with anatomical atlases), and in computer 
vision (target localization, automatic quality control). 

Many techniques of the image registration have been 
proposed previously. Medha V. Wyawahare et al. [1] 
proposed Image registration as a vital problem in medical 
imaging with applications in clinical diagnosis (Diagnosis 
of cardiac, retinal, pelvic, renal, abdomen, liver, tissue etc 
disorders). Ashok Veeraraghavan et al. [2] presented an 
approach for comparing two chains of deforming shapes by 
both parametric models like the autoregressive model and 



autoregressive moving average model and nonparametric 
methods based on Dynamic Time-Warping. Kendall's 
definition of shape is used for feature extraction. E.D. Castro 
and C. Morandi [3] placed steps for the new technique, based 
on the Fourier-Mellin transform. L. G. Brown [4] stated it as 
an intermediate step. J. P. Lewis [5] described it as cross 
correlation efficiently executed in the transform domain. S. 
Zokai and GWolberg [6] gave an idea about registration of 
multi resolution and multisensory images as a challenging 
problem in the research area of remote sensing and proposed 
new method for automatic affine image registration based 
on local descriptors, called automatic image registration by 
local descriptors (AIRLD). D. Lowe [11] presented a 
technique for extracting distinctive invariant features from 
images and achieves reliable matching among different views 
of an object or scene. J. B. Antoine Maintz et al. [12] 
presented a review of recent publications concerning medical 
image registration techniques. In probability and information 
theory, Mutual Information a dimensionless quantity with 
units of bits (sometimes known as transinformation) between 
two discrete random variables is termed as the amount of 
information shared between the two random variables. It can 
reduce the uncertainty. High Mutual Information (MI) 
indicates a large reduction in uncertainty; low MI indicates 
a small reduction; and zero MI between two random variables 
means they are independent. 

^-S^feSfe)' CD 

X and Y - Two discrete random variables. 
p(x,y) - Joint probability distribution function of 
X and Y. 

pl(x) and p2(y) - Marginal probability 
distribution functions of X and Y respectively. 

II. IMAGE REGISTRATION USING APT 

A. Localization 

The translation parameter between the two images was 
found. Fourier phase correlation is used to fix the translation 
before calculating the log-polar matching in the frequency 
domain. A study in [14] showed the aliasing problem by using 
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the phase magnitude of Fast Fourier Transform (FFT) to pick 
up the translation (of the centre point). 




Fig 1 . Feature point extraction and localization procedure for 
Transaxial Tl MR JPEG image of size 718 kB with resolution 1044 
x2219 pixels. 

First, feature points are extracted from both the 
model and target image of a single modality, either CT or 

MRI. P" = {p* Pj w } and F T = fof p& } are sets 

of feature points in the model and target image, respectively. 
The superscription and M and T denote the model and target 
image, and n M and n T are the number of feature points in the 
model and target image, respectively. To crop the model 
image out of the image, automatically one of the feature 
points is selected in the model image p£* e as a centre point 
for APT and chooses the size of the radius for the 
transformation R max that adequately covers the model image 
to compute APT and E M the projections of APT, R M and 
E M . On the left is the model image (Brain Image with the 
size 1044 x 2219 pixels) with 166 feature points identified. 
The model image is created by selecting one of the feature 
points in the model image pfas the origin, and is cropped to 
a circular image patch that covers the area R desired to be 

L max 

registered to the target image. On the right is the target image 
(Brain Image that is scaled by 0.7 and rotated by 30 degrees, 
the size of the image is 730x1553 pixels) with 180 feature 
points identified. To obtain the translation between the model 
and target image, the corresponding feature point p£J is 
located in the target image. As shown in Fig 1, the search 
space is limited to only a set of feature points in the target 
image, which is approximately 0.2% of the total number of 
pixels in the target image in this example. Hence, the search 
space is much smaller than that of the exhaustive search 
strategy. 

B. APT Matching 

After extracting feature points in the model image 
and selecting one feature point pj* as the centre point for 
computing the projections R M and E M using the APT 
approach, the next step is to find the related feature point 
p? r in the target image and obtain both the scale and rotation 
parameters between the two images of the same modalities. 
Given a set of feature points in the target image and the radius 
size of the APT transform R to be identical as for 

max v 

computing the projections R M and E M , each feature point is 
applied in the target image as a centre point for computing 
APT creating the set of candidate 

projections^ = [R[ f fl£ r }and e T = {&[ Gl T }. By 

matching and with each member in the sets of projections 
R M and E M , respectively, the translation, scale, and rotation 
parameters are obtained concurrently. The matching results 



have three dimensions: the scale, rotation, and distance 
coefficient. The translation parameter is the offset involving 
the locations of to the feature point in the target image that 
defers the lowest distance coefficient. 

For a given feature point in the target image, I M S /: , 
R- and O . are denoted as the model image, the image area 
that is cropped from the target image to be matched with 
the model image, projection on the radius coordinate, and 
angular coordinate, respectively. Then the scale change in 
the Cartesian coordinates reflects as variable- scalene of the 
projection, and rotation change reflects as shifting in the 
projection. Therefore, the new matching mechanisms are 
put forth to obtain both the scale and rotation parameters. 
The mathematical expression of image /lie, scaled and 
rotated version of image I^Jls 

l\ Of. y) = I M (a z xcos8 z + aysin& 2 . -a z xsin& 2 -f a g ycos& 3 ) (2) 

In the polar domain, (2) can be expressed as 

/ J o. ey = t* t ia. z *r*& —e x y 0) 
Where a^ and Q x are the scale and rotation parameters, 
respectively. As shown in (2) and (3), the scale parameter 
a z , in the Cartesian coordinates appears as scaling with the 
same value in the radius direction of the polar transformed 
image. For the case of adaptive sampling in APT, the 
outcome of scaling is seen in the Cartesian to APT mapping 

remains similar to that of the polar transform mapping. 
Except for APT, the number of samples in the angular 
direction can vary (adaptive). To calculate the ID projection 
in the radius direction R(i), multiplication variable U. is 
initiated to compensate the effect of adaptive sampling in 
the angular direction. As a result, the ID projection R of 
APT is an approximation to that of the polar transformed 



image. 

l£G>= 0*1^11? CM) (4) 
and fi H G)=n i 2£lP M G.j) (5) 
is obtained. R T _ is calculated in terms of R M as follows: 

fiJW^n^i^n^c^M) (6) 

JTJCO = R M (a,0 (7) 



Thus, variable-scale a z is a global and uniform scale 
for the ID projection R curve. As shown in Fig. 2(a) and (c), 
the projections E of the scaled images in the Cartesian is a 
little altered compared with that of the original image because 
the areas enclosed during the APT transforms are unlike result 
of scaling. Hence, to obtain the shifting parameter involving 
the two projections E M and E T and, variable-scale parameter 

a z involving the two projections R M and R T _ must be obtained 
first. 

B.A. Find Scale Parameter: 

Depending on prerequisites of the application in terms 
of the computational cost, accuracy, image types and 
environments, two effective algorithms are initiated. 
Algorithm I: Fourier Method: The Fourier transform of the 
scaled signal x (t) can be expressed as 

*Ct)~*ki) (8) 
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Similar to polar transform, the projection calculated 
from image /; that is the scaled version of image I M with the 
scale parameter a z will also be scaled version of the projection 
with the same scale parameter a z . From (8) and (9), given 
Fourier transform of the projection R(.) as A(.), the scale 
parametei ■ retained by /77 . 

This Fourier technique executes well only when the scale 
parameter is small. Large scale parameter defers the aliasing 
problem. Hence, this algorithm is apt for applications that 
involve fast and accurate image registration while the scale 
parameter involving the two images is expected to be small, 
such as medical image registration, or scene change 
detection. 

Algorithm 2: Logarithm Method: The second algorithm 
applies the scale invariant property of the logarithm function. 
First, the logarithm function is applied to the projection R 
and the output is then quantized to maintain the original 
dimension of the projection. The mathematical expression 
of the operation is as follows: 

LR(k) =R(n r ^ L )\k =1 n r (12) 

The parameter LR (11) shows the logarithmic of 
the projection. Given image r_ a scaled and rotated version 
of image I M , the scale parameter a z involving the two images 
would shows translation in the logarithm domain 

£LRl(p) = LR M (p) + logo, (13) 
To find the displacement d, where d=log a z , such 
that.£i?J (p) = LR M (p — d~) , the correlation function can be 
assessed between the two logarithmic of projections 

d = argmaxC(£R M ,£Rl~) (14) 
This algorithm defers high accuracy and fast matching in 
the general condition. 

B.B. Find Rotation Parameter: 

Next step is to find the rotation parameter e z . For the same 
sampling radius, R max , the larger the image (a z »l), the 
smaller the area of the scene is covered in the sampling 
procedure. As a result, the magnitude of the projections E 
between the two images could be slightly altered and detail 
is shown in Fig. 2(a) and (c). As a result, to precisely acquire 
the rotation parameter e z , (9) must be adapted according to 
the scale parameter a z : examine (14) and (15). Both 
projections are then resampled to be equal in length. The 
rotation parameter can be found by evaluating the correlation 
function£(0 >f , J) 



target image can be worked out from the Euclidean distance 
between the projections S M mnd Oj 



2nd 

"6 



(16) 



B. C Find Distance Coefficient: 

The final important component resulting from APT 
matching is distance coefficient, denoted as a z , which shows 
how large the variation is that involves two images. The 
distance coefficient a z between the model image I M and the 



(17) 



For registration, a feature point is searched in the 
target image that corresponds with the feature point in the 
model image. Such feature point is the point that yields the 

lowest distance coefficient. The projections R M and R T _ are 
not used for computing the distance coefficient a z as the 
dimension of projections that can be utilized in the 
computation varies depending on the scaling parameter a z . 
So, only projections § M and § - are used in the computation. 

C Image Comparison 

Using advantages of APT, quick and simple comparison 
method is put forth that can effectively and automatically 
position the altered area or area where the registered image 
pair differs without further computational cost. Since the 
scale difference in the Cartesian coordinates between the 
model image and target image yields different area coverage 
in the transformation process, the dimension of the projection 
R for comparing the two images needs to be altered 
accordingly. Given the scale and rotation parameters between 
the model and target image as a k r and e k , respectively, the 
altered area is found as a set of pixels (X, Y) in the Cartesian 
coordinates, where X=xl , . . . .,x , Y=y 1 , . . . ,y , and parameter 
is number of pixels in the altered area as follows: 

*«> HKr»-^{^LV^%Vi < I8 > 
£ e Q) = llSfW - S-0 - d)\\.f°rj = i to * 9 {19) 




HgJL Example of the proposed image registration approach (a) model 
image where the model image is the area inside the circular line, (b) 
two candidates are chosen as examples using two feature points A and 
B of the target image, (c) APT of the model image, (d) and (e) APTs of 
the target image where the centre points of the transformation are 
feature points A and B : respectively. 

T 

The location of p". . is denoted in the Cartesian coordinates 

Sir 

as (c x ,c y ). ForVff,. (j) , Vs e (j) > e, set of image pixels is 
computed in the altered area (X,Y) as follows: 

CYT) = {t lJ y 1 l....,Cx ff y t )} (20) 



(21) 
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The altered area can be recognized without extra 
computational cost. 



;?■■■ 



otherwise 



^[%IPi Ucrt (ti)J + flJ/J? ff.crf (£)Y «k#ru- !s . 

Most registration algorithms, involving the 2D 
geometric transforms, additional computations are essential 
as post processes such as image alignment; transform 
both.Images to the same coordinates prior to the comparison, 
and image normalization. For the proposed approach, as 
both the model and target image are already transformed 
to the projections of APT domain, image comparison can be 
performed directly and simultaneously while obtaining scale 
and rotation parameters. 

III.IMAGE REGISTRATION USING MAPT 

A multiscale transform of a signal U may be defined as 
I/fo, a) = J ttfr)^fr,fc a) dx t (24) 
Where a > is the scale parameter, b is the time (or position) 
parameter. 

Table I. 

Results for APT and MAPT for Brain image 
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In particular, analyzing function of Short-Time Fourier 
Transforms (STFT) using a Gaussian window (also known 
as Gabor transform). 



F^rfo b, a) = 9 * 9^- b \ 



(25) 



where j = V— 1- The analyzing function in wavelet 
transform using a Gabor wavelet (also called a 
Morletwavelet) is considered. 



(26) 



Where u Q is the basic frequency of Gwr The Gabor and 
Gabor-wavelet transforms (GWT) are implemented to do 
local frequency analysis as they are well localized in the 
time and frequency domains. The Gabor and GWT analyze 
signals at different frequencies and scales in a dependent 
way, i.e., the frequency content range depends on working 
scale. This is undesirable for image analysis as the coding 
regions are characterized by a fixed periodicity that may 
occur at different scales. Neither of these transforms is 
suitable for such task. 



b,b,a) = b~ B i«ot*-!0 (27) 



Therefore, Modified Gabor Wavelet Transform (MGWT) 
is defined as a function of b and a as 



(28) 
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IV. CONCLUSION 

A series of experiments was performed using 
medical images different sizes. Registration of monomodal 
brain images using MRI- Tl, T2 and CT scan images are 
studied and examined with regard to accuracy, robustness. 
For the algorithm evaluation, 41 such sets of CT image pairs 
were used. The results using MGWT in MAPT showed that 
scaled and rotated images registered well than the translated 
images. Further step is to implement a method that improves 
mutual information for translated images also. 
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